The 21-cm Signature of the First Stars 
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ABSTRACT 

We predict the 21-cm signature of the first metal-free stars. The soft X-rays emitted by these 
stars penetrate the atomic medium around their host halos, generating Lya photons that couple 
the spin and kinetic temperatures. These creates a region we call the Lya sphere, visible in 21-cm 
against the CMB, which is much larger than the HII region produced by the same star. The 
spin and kinetic temperatures are strongly coupled before the X-rays can substantially heat the 
medium, implying that a 21-cm absorption signal from the adiabatically cooled gas in Hubble 
expansion around the star is expected when the medium has not been heated previously. A 
central region of emission from the gas heated by the soft X-rays is also present although with 
a weaker signal than the absorption. The Lya sphere is a universal signature that should be 
observed around any first star illuminating its vicinity for the first time. The 21-cm radial profile 
of the Lya sphere can be calculated as a function of the luminosity, spectrum and age of the 
star. For a star of a few hundred M Q and zero metallicity (as expected for the first stars), the 
physical radius of the Lya sphere can reach tens of kiloparsecs. The first metal-free stars should 
be strongly clustered because of high cosmic biasing; this implies that the regions producing a 
21-cm absorption signal may contain more than one star and will generally be irregular and not 
spherical, because of the complex distribution of the gas. We discuss the feasiblity of detecting 
these Lya spheres, which would be present at redshifts z ~ 30 in the cold dark matter model. 
Their observation would represent a direct proof of the detection of a first star. 
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The formation of the first stars marks the end 
of the cosmic dark age and the beginning of the 
reionization era. The cold dark matter model 
with a cosmological constant (ACDM) of struc- 
ture formation predicts that the first stars were 
formed in dark matter halos of ~ lO 6 A/0, which 
were sufficiently massive to allow the halo gas to 
collapse and cool through emission in molecular 



hyd rogen roto- vibrational lines (for reviews see , 
e.g..lBarkana fc Loebll200lHMiralda-Escuddl20o3: 
ICiardi fe Ferraraf 2005 ). Numerical simulations 
and analytic arguments suggest that in the ab- 
sence of heavy elements, the gas in the halo center 
will cool hydrostatically and collapse into a central 
object when the Jeans mass is ~ 10 3 M Q , proba- 
bly resulting in th e formation of one or a few stars 
with ~ lOOM^ rtBromm. CoppL fe Larsonl | l99S 

Ja|2 



lAbel. Brvan. Normarj200dlOrnukai fe Pallaj2 003) . 
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20011 ) . At the end of their lives the supernova 



Stars of this mass live for only 3 x 10 6 yr ( Bromm. Kudri tzki HirEt^l ioodt. Pritchard fc Furlanetto 20061 ) . At 

the same time, the ionizing photons from the same 
first sources start heating the neutral medium, 
due mainly to soft X-ray s which penetrate deeper 
into t he atomic r egions ( Chen fe Miralda-Escudel 
2004 ISethil 120051 ). The 21-cm fluctuations are 



explosions enrich the intergalactic medium with 
heavy elements, which accelerate the cooling pro- 
cess and allow less massive stars to form. The 
reionization of the medium advances as more 
sources of ionizing radiation (stars or black holes 
emitting light as quasars) are formed in increas- 
ingly massive halos. 

A very powerful observational probe of the end 
of the dark age is 21-cm tomography, which pro- 
vides information on the state of the atomic inter- 
galactic medium (IGM ) fMadau. Meiksin. fc Reed 
19971: iTozzi et al.ll2000L fo r a recent review see e.g. 
Furlanetto. Oh fe Brigga 120061 ). Hydrogen atoms 
interact with 21-cm photons by radiative tran- 
sitions between the hyperfine split ground state. 
The hydrogen may be seen in emission or absorp- 
tion against the Cosmic Microwave Background 
(CMB) depending on whether its spin tempera- 
ture is above or below the CMB temperature. The 
spin temperature is a weighted average of the gas 
kinetic temperature and the CMB temperature: 
it is driven toward the CMB temperature by the 
21-cm transitions, and towards the kinetic temper- 
ature by atomic collisio ns and by resonan t scat- 
terin g of Lya photons ( Wouthuvsenl Il952t iFieldl 
1959). During the dark age (before the appear- 



ance of sources of light such as stars or accreting 
black holes), the kinetic temperature drops be- 
low the CMB temperature because of adiabatic 
cooling, and collisional coupling of the spin and 
kinetic temperatures can produce an absorption 
signal ( Loeb fc Zaldarriagal2004l : lBarkana fc Loeb 
2005bt iNaoz fc Barkanal 120051 ) . Collisional cou- 
pling becomes ineffective below a redshift ~ 40, 
except in high density regions such as colla psed 
minihalos (jlliev et alll2002t Kuhlen et alll2006l ). At 
lower redshift, the gas may again appear in absorp- 
tion when illuminated with the Lya photons from 
the first sources, which also couple the spin and 
kinetic temperatu res through the Wouthuysen- 
Ficld mechanism ( Chen fc Miralda-Escu dc 2004; 



1 This picture of dark age would be drastically changed 
if non-conventional source of light, e.g. decaying 
dark matter, were present. For discussions on the 
state of IGM and 21cm sig n al in such cases, see, e.g., 
IChen fe Kamionkowskilll2004l):lFurlanetto. Oh fc Pierapo'lil 



12006ft; 



IValdes et al al 



Ripamonti. Mapelli. Ferraral J2007f) 7 

ll2007t) . 



then due to variations o f neutral density, Lya flu x 
and spin temperature ( Barkana fc Loeb! 2005al ). 
As reionization proceeds, X-rays heat the ki- 
netic temperature to values much higher than the 
CMB temperature, and at the same time a Lya 
photon background is produced which strongly 
couples the spin and kinetic temperatures. The 
21-cm fluctuations are then essentially propor- 
tional to the density of neutral hydrogen, depen- 
dent on the IGM density and the ionization state 



( Madau. Meiksin. fc Reesl Il997 



Ciardi fc Madau 
Furlanetto. Sokasian. fc Hernquistl 



2003; 
Kohler et alll2005h . 



2004 



In the present work, we study the 21-cm sig- 
nature of the first stars, when the region around 
a star is illuminated by Lya photons for the first 
time after the recombination epoch. Before a Lya 
photon background is created, the spin tempera- 
ture is close to that of the CMB temperature at 
most spatial locations. However, near a first star, 
the Lya photon flux is strong enough to change 
the spin temperature, rendering the medium ob- 
servable with 21-cm tomography. This region is 
substantially greater than the HII region gener- 
ated by the star and the virial radius of the star- 
forming halo, as we shall see below. We shall refer 
to such a region as a Lya sphere. The gas in the 
inner part of a Lya sphere is heated up quickly by 
the radiation of the star (which, owing to the sur- 
face temperature T e ff ~ 10 5 K of massive, metal- 
free stars, emits substantial soft X-rays during the 
main-sequence lifetime) to a temperature above 
Tcmb- Therefore, this inner region produces in- 
stead a core of 21-cm emission. The heating rate 
is much lower at the outer part of the Lya sphere. 
Depending on the spectrum and the age of the 
star, the gas in the outer part of the Lya sphere 
should appear as a halo of 21-cm absorption. This 
outer halo of absorption is not present if the re- 
gion has been illuminated before by a sufficiently 
strong background of soft X-rays, and therefore 
heated by other more distant sources prior to the 
appearance of the star. Therefore, the detection 
of the absorption halo provides a unique signature 
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of a first star. 

As this work was being completed, we be- 
came aware of the work by Cen (2006) along 
similar lines. Our conclusions differ from those 
of Cen in two respects: first, we find that Lya 
spheres of absorption would be produced only 
around the very first sources of light in the uni- 
verse, which should be metal- f ree stars (see, e.g. 
Abel. Bryan, fc Norman] I2002L Bromm fc Larson! 
20041 ). At later times, the medium should be 
rapidly heated by X-rays and a global Lya back- 
ground should be produced, eliminating the Lya 
spheres. Second, the most important source of 
Lya photons that couple the spin and kinetic tem- 
peratures around a first star are not the continuum 
photons from the stellar photosphere, but the sec- 
ondary Lya photons due to impact excitations by 
electrons produced by the X-rays themselves. As 
we shall see in this paper, this fact is essential for 
improving the observability of these Lya spheres 
in 21-cm absorption, by increasing their contrast 
relative to the background. 

In §2, we derive the properties of an individ- 
ual Lya sphere. In §3 we investigate the forma- 
tion history of the first stars. The growth of the 
Lya and soft X-ray backgrounds, which affect the 
global evolution of the kinetic and spin tempera- 
ture of the IGM, is considered in §4. In §5, we 
consider the observability of these Lya spheres. 
We summarize our results in §6. In the following, 
we have adopted the WMAP three year ACDM 
model, with h = 0.73, O m = 0.23, fl b = 0.047, 
n s = 0.99, cr 8 = 0.77. 

2. Properties of Lya spheres 
2.1. Basic concepts and notation 

The form of the Lya sphere produced around a 
first star, detectable by means of the 21-cm emis- 
sion and absorption, is determined by the intensity 
of Lya photons and soft X-rays produced over the 
history of the stellar lifetime. As a result of the ir- 
radiation of the neutral medium by stellar ultravi- 
olet light, Lya photons are produced by two mech- 
anisms: they may have been emitted directly from 
the stellar photosphere at a frequency in the range 
between Lya and Ly/3 and then redshifted to Lya 
, or they may result from ionizing photons which 
can ionize hydrogen atoms that subsequently re- 
combine and also excite other hydrogen atoms 



with the energeti c electrons created in the ioniza- 
tions. Following fChen fe Miralda-Escudl (|2004l ). 
we designate the first type continuum photons, 
since they enter the Lya line from the continuum 
spectrum on the blue wing of the line, and the sec- 
ond type injected photons, since they are injected 
close to the line center by an atom that has recom- 
bined or has been excited by a secondary electron, 
and has then reached the 2p state. We note that 
injected photons may also originate from photons 
emitted between the Ly7 line and the Lyman limit 
by the star, which will be absorbed by a hydrogen 
as they reach a Lyman resonant line, and may 
then produce a Lya photon if the hydrogen atom 
cascades down to 2p. However, the majority of 
these photons result in a cascade that ends in the 
2s state instead, so these typically make a small 
contribution to the creation of injected Lya pho- 
tons compared to the ionizing photons (see Hirata 
2006; Pritchard & Furlanetto 2006). 

A fixed point in the intergalactic medium can 
be reached by continuum photons from stars out to 
a maximum redshift of z c ^ max — (8/9)/ (3/4) — 1 = 
5/27. Injected photons from ionization may come 
from a higher redshift, but they need to arise from 
X-rays which can penetrate over a large distance 
in the atomic medium. Every X-ray can produce 
several secondary ionizations and collisional exci- 
tations by the electrons that are produced, and 
these will result in the production of a Lya photon 
when the hydrogen atom reaches the 2p (rather 
than the 2s) state. These Lya photons start scat- 
tering after they are produced, and they are con- 
fined to a very small region of the universe owing 
to the large cross section for Lya scattering in the 
neutral medium. The flux of Lya photons at any 
point in the atomic intergalactic medium is there- 
fore related to the local flux of X-rays and con- 
tinuum photons between the Lyman series lines 
coming from distant stars, although some spatial 
diffusion of the Lya photons may take place in 
the last fe w scatterings, when th e photon exits the 
Lya line (|Loeb fc Rvbickilll999l) . 

The brightness temperature of 21-cm photons 
is given by 



ST h = 40mK 



0.3 



0.03 ) \n m0 

. Phi T s — Tcmb 



1/2 



1/2 



PH 



(1) 
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The spin temperature T a is a weighted average of 
the CMB temperature and the gas kinetic temper- 
ature, and can be approximately calculated as (see 
Hirata 2006 for a more exact treatment) 



T s = 



Tcmb + {y a + Vc)T k 



(2) 



where the coupling coefficients y a an d y r are pro- 
portional to the Lya scattering rate (jWouthuvsen 
1952HFieldlll959h and spin-changing collision rate, 
respectively: 



PiqT* 
A w T k • 



A w T k ' 



(3) 



Here, T* — 0.068K is the temperature correspond- 
ing to the energy of the hyperfine structure of hy- 
drogen, A w — 2.87 x 10~ 15 s _1 is the Einstein co- 
efficient of the hyperfine structure levels, and Cio 
is t he collisional de-exc i tation rate (tabulated first 
bv lAllison fc Dalgarnol (Il969h. and more recently 
bv IZygelmanl (|2005l) ; IFurlanetto fc Furlanetto 
(|2007al lbl): in t his paper we use the tabulated 



values given bv IFurlanetto. Oh fc Briggsl ( 20061 )1. 
The spin-flip transition rate Pio is given by 



4 S C J C + 

PlO = 7^ Ht GP f 

Li Jo 



(4) 



where H is the Hubble rate at that redshift, tqp 
is the Gunn-Peterson optical depth, J c and Ji are 
the intensities of continuum and injected Lya pho- 
tons, and S c and Si are numerical factors of order 
unity tha t depend on the gas temperatur e and on 
top ( see IChen fc Miralda-Escudl (|2004l ); iHirata 
( 20061) Fll. In the high temperature limit, S c ,i — > 1. 



For simplicity, in the present paper we shall as- 
sume S c = Si, an d use the value as given by th e 



fitting formula of IFurlantto fc Pritchard (2006). 
The fiducial intensity Jo corresponds to a photon 
density of one photon per hydrogen atom per log 
frequency, 

Jo = -. , (5) 



2 The numerical values of the hea ting; and collisional rate 
coeffic ients Ii, I c , Si,S c given in IChen fc Miralda- Escudel 
are all too large by a factor of t/tt, because 
of an incorrect normalization used in the original code 
(this does not affect the ba sic conclusions reached in 
IChen fc Miralda-Escudl 1120041) 1. We thank Chris Hirata 
for pointing out this error to us. Also, loss of Lya photons 
by the two photon decay process could further reduce the 
effective collision rate slightly, sec Hirata (2006) for further 
information. 



where nn is the number density of hydrogen 
and v a is the frequency of Lya . In ordi- 
nary units, Johv a — 2.69 x I0~ 24 x (I + 
z) 3 (fi 6 ^/0.02)ergcm- 2 s" 1 sr- 1 Hz _1 . At z ~ 
20 - 30, y a ~ 1 when J/J - 0.02. 

2.2. Lya sphere model 

The Lya sphere of a first star is the region 
where the spin temperature is substantially modi- 
fied by the Lya photon flux arriving from the star, 
producing an observable signature in the 21-cm 
line. Stars that form later emit their ultraviolet 
light into a region that has already been previ- 
ously illuminated by a background of light from 
more distant stars which formed earlier, and the 
contrast of the Lya sphere from an individual star 
will of course rapidly diminish as the background 
intensity increases. In this paper, we designate as 
"first star" any star that emits the first light into 
the atomic medium that surrounds it. Many more 
metal-free stars form later from pristine material 
(i.e., not contaminated by any stellar ejecta), but 
which has been illuminated by stellar light; these 
stars create a Lya sphere in the presence of a light 
background, with lower contrast. Note that the 
region around a star is also modified by an in- 
creased kinetic temperature due to ionization by 
soft X-rays. 

In this p aper we assume that all the first stars 
are massive ( Bromm fc Larsonl 2004), with masses 
ranging from 25 iW© to 800 Mq. The main se- 
quence lifetime of a massive, metal-free star is 
about three million years, and is generally much 
longer than the Kelvin- Helmholz time to reach the 
main-sequence. We treat the stars as static dur- 
ing their lifetime, with constant luminosity. As 
we shall see below, due to the r~ 2 decrease of the 
photon flux, the size of Lya sphere is typically 
limited to a few tens of kpc, so the corresponding 
light propagation time is also negligible compared 
with the stellar lifetime. The scale of this Lya 
sphere is much greater than the halo harboring 
the star, so we treat the density of the gas as uni- 
form and equal to the cosmic mean (note that this 
may be inaccurate in many cases because of the 
highly biased large-scale distribution of the first 
stars) . The flux of the stellar continuum per unit 
frequency is 

L(v a ) 



Jr. = 



(4tt) 



(6) 
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where r p is the physical distance to the star, and 
L(v) is the stellar luminosity per unit frequency. 
We model the stellar radiation as blackbody with 
surface temperature and total lum inosity given by 
I Bromm, Kudritzki fe Loebll200lh 



T = 1.1 x 10' 



M 



( M 



V 100M F 



0.025 



Ltotal — 10 4 '°— — 

Mq 



K (7) 
(8) 



2.3. 



Ionization and heating of the Lya 
sphere 



We now proceed to solve for the ionization frac- 
tion, kinetic temperature, spin temperature, and 
the 21-cm emission temperature around the star. 
We include recombinations in our calculation, but 
as the recombination time at the average density of 
the universe is a few times longer than the lifetime 
of the star and its associated Lya sphere, it does 
not have a big impact on our result. Note that we 
are also assuming that the high-density halo gas 
around the star (which does not have time to be 
ejected from the halo during the stellar lifetime af- 
ter it is ionized and heated) does not absorb a large 
fraction of the st ellar ionizing photons emitted b y 
the central star (jWhalen. Abel, fc Norman 12004 . 
We also do not consider ionization and heating by 
any X-ray background that may have been pro- 
duced by sources other than the star. The neutral 
fraction is determined by 



d.X 



HI 



dt 



dt 



d.r 



Hell 



dt 



/>oo 

-xm / dv F(v)am{v) + a^xmxenn 

J V H 

(9) 

-x Hc i / d^F(t/)(THci(^)+ai eI a:Hoiia; e «Hc 

(10) 

XHeii [ dvF(v)aK c i\{v) 



dt 

+a^ cII X e (l - XHoI - XRcIlWc (11) 



Here, F(y) is the flux of ionizing photons per unit 
frequency 7jh and riHe are the total number den- 
sity of H and He atoms, and Xi , (y) , and Vi are 
the fractional abundance, the cross section and 
the ionization threshold frequency of each atom 
or ion of type i, respectively. The case B recombi- 
nation rates for hydrogen is calculated by us- 
ing the fitting formula given in lHummer fc Seaton 



(|1963l ), the ag , ctg are calculated by inter- 



po lating/ extrapolating the ta b ulated values given 



Storey fc HummerT(|l995h : iHummer fc Storevl 



( 19981 ). In the helium ionization, we have ne- 
glected the contribution of direct double ioniza- 
tion, as this rate is small compared with ionization 
in two stages. We fit the cross s ection s using the 
functions given in IVerner et all (|l996h . The flux 
at physical distance r p is given by 



F{y) 



(12) 



where 



nB_o m (v) I dr'x m (r') 
Jo 

+«h c cthci(^) / dr'x Ue i(r') 



+«HeO-HeIl(^) / ^'^HellC^') (13) 
JO 



2 X 10 4 , X H eII 



We assume xhi = a;Hoi 
everywhere initially (to include the primordial 
ionization). We then evolve the ion fractions on 
a lattice of 200 points (logarithmically distributed 
from r = O.Olkpc to r — lOOkpc) with the above 
differential equations. The solution of the ordinary 
differential equati on is obtained with a predicto r- 
corrector scheme ( Engeln-Mullges fc Uhlid [T996). 

The internal energy of the gas is given by 



u = 7; ["h(2 - xm) + «Hc(3 - 2x Hc i 



EHell)] kT 

(14) 



The volume heating rate r = ^ is given by 



i=HI,HeI,HeII 



dvh{y—Vi)r}{x, v)F(v)a l {v), 

where the fraction of photo-electron energy con- 
verted to heat J](x, v) is a function of energy 
and fr ee electron fraction. IShull fc van Steenberd 
(|1981l ) calculated this function using a Monte 
Carlo simulation, and found that it can be fitted 



a.s 



rj = 0.9971(1 - (1 



„0.2663\1.3163 



^1.31B3^ 



(16) 



Note that this fit breaks down at very low x, but it 
is a good approximation for x > 10" 4 . Also, this 
fitting function was obtained for non-primordial 
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helium abundance, so the real result may differ 
slightly. 

At large r p , where the ionized fraction in the 
Lya sphere is low, cooling is negligible and we 
simply calculate the temperature using the above 
expressions for calculating the internal energy and 
the heating rate, assuming that there is no cooling. 
The effect of atomic cooling sets in rather abruptly 
once the temperature reaches T ~ 10 4 K, and we 
assume that the temperature stays fixed at this 
value after it is reached. This does not affect our 
result on 21-cm because at these high tempera- 
tures the emission depends only on the hydrogen 
column density. 

2.4. Secondary Lya photons from ioniza- 
tion by X-rays 

Owing to the high temperature of the metal- 
free stars, the soft X-rays from the stellar photo- 
sphere result in a substantial ionization and heat- 
ing ahead of the thin shell where most of the ion- 
ization occurs. These high-energy electrons pro- 
duced by the X-rays that penetrate into the mostly 
atomic zone also result in the emission of addi- 
tional Lya photons: in the limit of low fractional 
ionization and high photon energy, a fraction i] a ~ 
40% of the X-ray energy is converted to photons 
through excitation by the high-energy electrons. 
Here we assume all of these are in Lya photons. 
In reality, the fraction may be slightly smaller. The 
intensity of these Lya photons is added to those 
coming from the stellar photosphere and can con- 
tribute to the coupling of the spin and kinetic tem- 
peratures. 

The Lya photons produced through X-ray ion- 
ization are "injected" photons, since they are in- 
troduced at the resonance line center. Spatial 
diffusion of these photons is negligible, because 
the width of the scattering line is much less than 
the Hubble expansion velocity of the Lya sphere 
at the radius where the emission of these pho- 
tons is important. To evaluate their intensity, we 
note that each X-ray will produce a number of 
Lya photons rj a v/v a , and that each Lya photon 
produced in this way will have its frequency red- 
shifted at a rate Hv a if the effect of scatterings 
are neglected, where H is the Hubble constant. 
Hence, the number density of these Lya photons 
per unit of frequency should be equal to the heat- 
ing rate per unit of volume (eq. 15) times the factor 



{jja/vi) / {hv a ) / {Hv a ). Their intensity is therefore 
equal to 

c Va T 

In terms of the fiducial intensity in eq. (5), we can 
rewrite this as 



Ji 



J r]hv a n H H 



(18) 



Noting that the quantity hv a nnH is roughly the 
heating rate required to heat the gas to a tem- 
perature hv a /kB ~ 10 5 K over a Hubble time, 
we can easily see a relation between the temper- 
ature to which the gas is heated by X-rays and 
the Lya intensity that can be produced from the 
same X-rays. For example, considering the gas 
temperature and Lya intensity at half the lifetime 
of the star, t = 1.5 x 10 6 yr ~ lO" 2 !?" 1 , at a 
distant point in the Lya sphere where the tem- 
perature has been heated due to X-ray absorption 
by AT = 10 K (still allowing for a large 21-cm 
absorption signal), and using r/ a /rj ~ 4 (valid for 
high neutral fractions and X-ray energies), an in- 
tensity Ji/Jo ~ 4 x 10 -2 can be generated, suffi- 
cient to obtain y a ~ 2 for the spin-kinetic temper- 
ature coupling. 

We note here that the scattering of Lya photons 
can also change the kinetic temperature of the gas. 
As shown by Chen & Miralda-Escude (2004), the 
heating rate due to the Lya scattering is negligible 
compared to that caused by X-rays, and we do not 
include it here. In fact, for the regime that we shall 
be interested here, the temperature of the atomic 
medium is high enough that the scattering of Lya 
photons actually causes cooling of the gas (see Fig. 
3 in Chen & Miralda-Escude 2004), although by a 
negligible amount. 

2.5. Results on the ionization and temper- 
ature profile evolution 

We plot the neutral and ionization fraction of 
the gas as a function of distance from the star in 
Fig. □ for a star of 2OOM at z = 20 and age 150 
Myr. As can be seen from the figure, outside the 
HII region, the ionization fraction drops, but there 
is still some residue ionization, which is produced 
by the penerating X-ray photons. At about 15 
kpc, the ionization fraction approaches the back- 
ground value. The spin and kinetic temperature 
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Fig. 1. — This shows the neutral and ionizatgion 
fraction of H for a 200M Q star at z = 20, with age 
1.5 Myr. 
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Fig. 2. — The kinetic and spin temperature (low 
panel) and 21-cm brightness temperature (upper 
panel) of a first star (same condition as above). 



profile, and the 21-cm brightness temperature is 
plotted in Fig. H 

The spin temperature is strongly coupled to the 
kinetic temperature throughout the zone of the 
strong temperature drop ahead of the ionization 
front. The injected Lya photons generated from 
the ionization by X-rays are far more intense than 
the continuum Lya photons originating from the 
stellar photosphere. The y a parameter (see eq.[3|), 
which indicates how strongly the spin and kinetic 
temperatures are linked by the Lya photons, be- 
comes of order unity just at the radius where the 
atomic medium starts being heated above its ini- 
tial temperature by the X-rays from the star (for 
the reason explained above in §2.4). At larger ra- 
dius, the spin temperature gradually rises above 
the kinetic temperature towards the CMB tem- 
perature as the Lya coupling becomes weaker. 
This produces a region of strong 21-cm absorption 
against the CMB, coming from the atomic medium 
that is starting to be ionized and heated, still far 
from the ionization front where most of the ion- 
ization occurs. The temperature in this region has 
not yet increased very much, allowing for a strong 
absorption signal, but the Lya photons generated 
by the ionizations are already strong enough to 
lower the spin temperature by a large factor below 
the CMB temperature. Note that at very large ra- 
dius, the spin temperature reaches a constant that 
is slightly below the CMB temperature because of 
collisional coupling, causing a constant level of ab- 
sorption on the CMB. 

Inside the HII region (the region that is mostly 
ionized), the kinetic temperature is assumed to 
stay constant. The spin temperature should also 
stay high, not only because of the Lya photons 
produced by the high-energy electrons generated 
in the ionization front, but also because of addi- 
tional Lya photons generated by recombinations 
in the HII region. We have not calculated the 
detailed shape of the spin temperature profile in 
the HII region (this would be affected by diffusion 
of Lya photons generated at the ionization front). 
However, this profile does not actually matter very 
much because, in the regime where the spin tem- 
perature is much higher than the CMB tempera- 
ture, the observed 21-cm brightness depends only 
on the column density of hydrogen per unit of ve- 
locity, and not its spin temperature. 

To see how these depends on the redshift, age, 
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Fig. 3. — Neutral fraction as a function of radius 
in the HII region at the indicated redshifts and 
times after the birth of the star. In each panel, 
the six curves are for a star of mass M — 25, 50, 
100, 200, 400, and 800 M , from bottom to top. 



and mass of the star, we made plots of Fig. [3] 
and Fig. 2] In Fig. [31 the top, middle, and bot- 
tom panels correspond to a time after the birth 
of the star of 0.5, 1.5, and 2.5 Myr. The result is 
shown at redshifts z = 30, 25, 20, 15 from left to 
right. In each panel we plot the neutral fraction 
for six stellar masses: 25, 50, 100, 200, 400, and 
800 Mq. The ionized regions grow with the age of 
the star; for fixed stellar mass and age, the HII re- 
gion is smaller at higher redshift because the gas is 
denser. The maximum physical radius reached is 
1 to 10 kpc. The kinetic and spin temperatures of 
the Lya sphere are plotted as a function of phys- 
ical distance in the lower panels of Fig. [4] at a 
time t = 1.5 Myr and the four redshifts z = 30, 
25, 20, and 15. The six curves in each panel are 
for stellar masses M = 25, 50, 100, 200, 400, and 
800 Mq. We have assumed that there is no global 
Lya background or heating prior to the formation 
of the star. The kinetic temperature is equal to 
our imposed floor (10 K) within a few kilopar- 
secs, due to heating by ionization. Beyond this 
radius the kinetic temperature drops rapidly as 
the atomic zone is penetrated, until it reaches the 



physical radius (kpc) 

Fig. 4. — Brightness and temperature profiles of 
Lya spheres, at a time t = 1.5 Myr and the indi- 
cated redshifts. Bottom panels: Tk (dash lines), 
Tcmb (dotted line), and T$ (solid lines). The 
six curves are for a star mass M = 25, 50, 100, 
200, 400, and 800 M Q . Top panel: 21-cm bright- 
ness temperature fluctuation on the CMB at time 
t = 1.5 Myr. 



constant value of the unheated atomic medium far 
from the star, which is below the CMB tempera- 
ture owing to adiabatic cooling. 

In Fig. [51 we plot the temperature profile if the 
Lya photons are produced only from the star, and 
not by X-ray excitation. As can be seen from the 
figure, the coupling is much weaker so the spin 
temperature is not strongly coupled to the kinetic 
temperature, and the absorption signal is much 
weaker. 

The Lya sphere evolves as the gas is heated and 
ionized by the star. In Fig.[6l the same profiles are 
plotted at redshift z = 20 and three different ages, 
t = 0.5, 1.5, and 2.5 Myr. The features of the 
Lya sphere shift with time to larger radius and 
become shallower as the region where ionization 
and heating are important grows. A cross section 
map of the Lya sphere for a star of 200M Q , age 
1.5 Myr at z — 20 with no global heating is shown 
in the left panel of Fig. [ITJ The central HII re- 
gion where ST), ~ 0, the emission region (red) and 
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Fig. 5. — The temperature profile without X-ray 
induced photons. 
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Fig. 6. — Brightness profile of Lya spheres at z = 
20, for the three indicated ages after the star birth. 



absorption region (blue) is clearly seen. In real ob- 
servations, the absorption signal against the CMB 
at large radius is the dominant signature of the 
Lya sphere formed around a first star: it is of 
much higher brightness temperature and occupies 
a much larger area than the emission region close 
to the star that has been heated by the ionization. 
This small emission region could only be detected 
with very high angular and frequency resolution. 
In the absence of sufficient frequency resolution, 
the projected region of absorption on top of the 
small emission region would greatly overwhelm the 
signal. 

3. Formation History 

We now consider the formation history of the 
first stars in order to estimate the abundance of 
the Lya spheres discussed in the previous section. 
We calculate the comoving density of collapsed 
dark matt er halos with the halo m odel of Press- 
Schechter ( Press fc Schechterl \1974. and also the 
Shcth-Tormen (jSheth fc Tormenl I1999T 1 prescrip- 
tion. The density of collapsed halos of mass m 
at redshift z can be written as 



N(m,z) = —sf{s)-^- 
m dm 



(19) 



where s = 5l(z)/a 2 (M ), a(M) is the variance of 
the density fluctuation smoothed on a mass scale 
M with a top hat filter, and the critical linear 
overdensity is S c (z) — 1.686/£>(z), where D(z) is 
the linear growth factor. The function /(s) is 



sf{s) 




-s'/2 



(20) 



Here, s' = as, with A = 1/2, a = l,p = for 
the Press-Schechter (PS) mass function, and A = 
0.322, a = 0.707, p = 0.3 for the Sheth-Tormen 
(ST) mass functions. 

The gas in a dark matter halo can form a first 
star provided that the primordial gas in the halo 
is able to cool by means of the small fraction 
of molecular hydrogen that can be made at this 
early epoch, faster than the rate at which the gas 
is dynamically shock-heated during mergers and 
accret ion (e.g., lYoshida et al. ( 20031 ). Reed et al 



(2005)). We assume for simplicity that star for- 
mation proceeds in any halo with T vir > 2000K, 
where T vir is the gas temperature after virializa- 
tion and in hydrostatic equilibrium. For neutral 
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gas, the mass of the halo when the virial temper- 
ature reaches this value is 



M mm (z) = 1.17x10 



, O h 2 
0.147 



-1/2 



(1 + *) 



-i -3/2 



21 



T 



2000K 



(21) 



The lifetime of one of these massive, first stars is 
i* S3 3Myr, so only halos that have formed within 
the last t* would possess an active star. We as- 
sume here that only one metal-free star forms per 
halo, since the ionizing radiation from these mas- 
sive stars is ab le to ionize and expel all the gas in 
the halo (e.g., IWhalen. Abel fc Normanl (|2004h ). 
and any future stars will already form from metal- 
enriched gas (presumably with a mass function 
closer to the present one). 

After its main sequence life time, a first star 
may explode as a supernova and destroy itself 
completely; or explode as a supernova and left be- 
hind a black hole or n eutron star; or collapse into 
a bla ck hole directly (jWooslev. Heger. fc Weaver 
l2002i ). In the absence of continued supply of X- 
ray and Lya photons, the Lya sphere dims as the 
scattering Lya photons diffuse out. However, su- 
pernovae remnant and accreting black hole may 
become new source of light, the Lya sphere could 
be maintained, or even grow to much larger size. 
We will investigate these interesting possibilities 
in subsequent work. 

For simplicity, we assume that one star is 
present in each halo that has just formed with 
mass M > M m i n (z). At high redshift, when the 
fraction of matter that has collapsed in halos in 
which gas can cool is small, the number of halos 
of mass m that have formed over a time t* can be 
approximated as 



dz 



M min (z) 



dN(m, z) 

dm , 22 

oz 



where the cosmic time at z a is t(z a ) — t{z) — £*. 
In the redshift range 15 < z < 40, this can be 
simplified as 



dm (23) 
oz 



n* S3 (1 + z)H(z)U 

lM min {z) 

At low redshift, the approximation obviously fails 
as halo destruction by mergers becomes impor- 
tant. At high redshift n(m, z) evolves too rapidly 
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Fig. 7. — Number density of star-harboring halos 
of mass m per comoving Mpc -3 per log to. The 
four curves are at z = 15 (solid, blue curve), z = 20 
(dotted green curve), z — 25 (short-dash yellow 
curve), and z = 30 (long-dash red curve). 



and the integration over redshift needs to be done 
more carefully. 

We plot the mass function of halos formed 
within a time i*, calculated with the Press- 
Schechter formalism in Figure [7] In Figure [5J 
the comoving densities of star forming halos (solid 
curves) and first stars (dashed curve; this is equal 
to the density of halos formed within the past 
three million years) are shown, both for the Press- 
Schechter and Shcth-Tormen mass functions. The 
ST mass function predicts a higher abundance for 
rare fluctuations, so at very high redshift the ST 
result is substantially greater than the PS result. 
At lower redshifts the ST function fits N-body sim- 
ulation results better, although it te nds to over- 



predi ct the number of halos at z > 10 (jReed et al 
2003). The accuracy of the mass functions at 
such high redshifts is still largely unknown, dif- 
ferent conclusio ns were reached in different nu- 
merical studies (Tjang-Condell fc Hernquistl 200ll: 
Jenkins et ail |200lt iReed et all 120031: llliev et al 1 



2006; Zahn_elaU2006; Wa rren et al . 2006; L ukic et al 



20071 ). and in fact may also depend on the d efini- 
tion of the halo and its mass in the simulation! Cohn fc White] 
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Fig. 8. — Comoving density of halos and first 
stars. Solid curves: total density of star-forming 
halos. Dashed curves: density of star-forming ha- 
los formed in the past 3 Myr. Results are shown 
for the Press-Schechter (PS) and Sheth-Tormen 
(ST) halo mass functions. 
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Fig. 9. — Bias factor of star-forming halos for the 
PS function (solid curve) and ST function (dashed 
curve). 



20071 ). 

The comoving density of first stars peaks 
around a redshift z ~ 15. Our calculation be- 
gins to fail at this redshift for several reasons: 
the merger destruction rate of the halos becomes 
comparable to the formation rate, and equation 
2"2l is no longer accurate. Physical feedback effects 
suppressing the formation of metal-free stars and 
new star formation in regions containing metals 
may become important even earlier than z — 15. 
We do not consider these uncertainties in this pa- 
per, where we merely want to illustrate how the 
Lya spheres of the first stars may interact, so our 
calculation is stopped at z = 15. 

Halos harboring the first stars are strongly clus- 
tered because of the flatness of the Cold Dark Mat- 
ter power spectrum at small scales, which implies 
a high bias of the first objects in large-scale regions 
of high density. The bias of halos of mass M at 
redshift z c an be approximated as (|Mo fc White! 
19961 l2002h 
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Fig. 10. — Correlation function of star-forming ha- 
los, at redshifts z = 15, 20, 25, and 30 (from bot- 
tom up), calculated with the Press-Schechter func- 
tion. Note that the correlation function increases 
with redshift, because of the increasing bias factor. 



(24) 
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Fig. 11. — Expected number of star- forming ha- 
los as a function of distance. The four curves 
are z = 15 (solid, blue curve), z — 20 (dotted 
green curve), z — 25 (short-dash yellow curve), 
and z — 30 (long-dash red curve). A physical dis- 
tance of 10 kpc (the approximate size of the Lya 
spheres) is marked at each redshift. 



for the PS function, and 
1 



1 



Sc(0) 



0.5s 



10. 4 



J0.6 



/0.84 



0.14 



(25) 



for the ST function. We calculate the bias corre- 
sponding to the minimal mass for halos harbor- 
ing stars at each redshift. The result is shown 
in Fig. [SJ The bias grows very rapidly with red- 
shift, faster than (1 + z), implying that (perhaps 
counter-intuitively) the first stars are increasingly 
clustered as the redshift increases, at a fixed co- 
moving scale. 

The correlation function of biased halos is given 

by 

£ hh (r,z)*ib 2 £ lin {r,z) , (26) 
where the linear correlation function is 



t(r,z) = 



d 3 k 



P Un {k,z)e 



■{ kr 



(27) 



We plot the correlation functions in Fig. [TDJ The 
increasing amplitude of the correlation function 
with redshift is again a result of increasing bias 
factor with redshift. 



The average number of neighboring Lya spheres 
within a given radius r is (neglecting higher order 
clustering terms) 

r R 

N(R)=n(z) / dV [1 + £(r)] , (28) 



where n(z) is the number density of star forming 
halos. Note that this calculation is not very ac- 
curate, as non-linear effect may be important at 
such small scales, it is only meant to give an idea 
on how the Lya spheres are correlated. The result 
is shown in Fig. [TT] The average number of halos 
within a sphere of physical radius r = 10 kpc is 
~ 0.04 at z = 30 and ~ 4 at z = 15. This number 
increases with time because of the growing num- 
ber of collapsed halos, despite the decrease in the 
correlation function. Even if several similar halos 
are found within the radius of the Lya sphere, the 
formation time of the two stars may not coincide, 
thus the observed density of Lya sphere is f*N, 
where the duty cycle factor /* ~ n*/rihalo < 1- 
Moreover, feedback effects (e.g., from molecular 
photodissociation) of the radiation of one star may 
adversely affect the formation of the other one. 

The large biasing factor of the star-harboring 
halos implies that the first stars which produced 
Lya spheres may have an irregular neighborhood, 
with more collapsed structures on the scale of the 
Lya sphere than one would find in a random lo- 
cation of the universe. These density fluctuations 
will likely make the "Lya spheres" have a struc- 
ture that is not actually very spherical, but irreg- 
ular. 

4. The impact of Lya and hard X-ray 
backgrounds 

4.1. The Lya background 

In §2, we have computed the kinetic and spin 
temperature profiles of the atomic medium around 
an isolated first star. The largest signal this 
"Lya sphere" produces in 21-cm is the absorp- 
tion against the CMB at large distances from the 
star. This absorption can be detected against the 
background if the spin temperature of the atomic 
medium is coupled to the kinetic temperature only 
close to the star, but not far from it because of 
the absence of Lya photons; each star would then 
produce an isolated Lya sphere. However, as the 
number of first stars increases, a Lya background 
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will be generated which can couple the spin and ki- 
netic temperatures everywhere, producing a global 
absorption signal against which the regions near 
one star would be harder to detect. In fact, owing 
to the presence of the foreground, it is not possible 
to make absolute measurements of the 21-cm sig- 
nal, but only its spatial variations and non-smooth 
variation in redshift. In Figures l4l and [T7l the plot- 
ted emission (absorption) temperature is not the 
absolute value, but the difference its value far from 
any Lya sphere. 

The Lya background originates from the con- 
tinuum photons emitted by the stellar photo- 
spheres. This is because the injected photons re- 
sulting from the high-energy electrons generated 
by soft X-ray photons are emitted at the Lya line, 
and are rapidly redshifted as they leave the Lya 
sphere where they have been created. On the other 
hand, hard X-ray photons may penetrate the Lya 
sphere and produce some injected photons very far 
from the sources, but the total energy released in 
hard X-ray photons by the first stars is less than 
that in UV photons, and only a very small fraction 
of the hard X-ray photon is absorbed in the IGM 
to produce injected Lya photons. Because the in- 
jected photons vastly dominate the intensity over 
the continuum photons from the central star, the 
background of Lya photons is much fainter than it 
would be if continuum photons were the dominant 
ones. This allows for a lasting contrast of individ- 
ual Lya spheres in comparison to the background, 
because the Lya background intensity stays low 
even when a very large number of first stars are 
seen from any random point in the universe. 

At the same time, we need to take into account 
the background of X-rays produced by the stars 
and their global heating effect, which will raise the 
kinetic and spin temperatures and hence decrease 
the strength of the 21-cm absorption. 

To compute the background of the continuum 
photons near the Lya line, we assume a blackbody 
spectrum for the first stars. Let E(y) be the num- 
ber of photons emitted over all time per baryon 
and per unit of frequency in one of the halos that 
collapse to form the first stars. In a halo that 
forms a single, massive metal-free first star, the 
spectrum should be a blackbody with temperature 
T ~ 10 5 K, with little dependence on the stel- 
lar mass (Schaerer 2002). As more massive halos 
collapse and continue to form stars from metal- 



polluted gas, different stellar populations should 
form and emit more photons. For simplicity, we 
shall assume here that E(v) remains constant, 
that is to say, the total number and spectrum of 
photons emitted per baryon in collapsed halos is 
fixed at frequencies greater than Lya (which are 
the only ones that affect our result). In view of 
our present ignorance on the type of stars that 
will form in halos of different mass as the processes 
of collapse of structure in the universe, reioniza- 
tion and metal pollution proceed, this assumption 
seems as good as any other. 

The fraction of mass bound in star-forming ha- 
los is 

f°° dN(m,z) , , . 

F(z) = / m \-^—!-dm, (29) 

Jm^z) dm 

where dN ^™' is the comoving density of halos of 
mass m, and M m i n (z) is the minimum halo mass 
to form the first stars at redshift z. The star for- 
mation time and the main-sequence lifetime of the 
first stars are short compared to the Hubble time, 
so the emissivity can be approximately written as 

e(u,z) = n H E(v)F . (30) 

where F is the time derivative of the fraction of 
mass bound in star-forming halos, F, and we have 
assumed a constant E(y). 

As mentioned above, the only photons that con- 
tribute to the Lya background are those emitted 
by the stars in the continuum to the blue of the 
Lya line. Among these, only photons with v < vp 
can be redshifted to v a as continuum photons, and 
all photons of greater frequency will be absorbed 
at the Lyman (3 or higher resonance Lyman series 
lines and be converted to injected Lya photons. 
Here we will calculate only the intensity of the 
continuum Lya background, and we will need to 
bear in mind that the total background will actu- 
ally be a little bit larger because of the injected 
photons from higher Lyman series lines (note that 
the spin changes caused by these higher Lyman 
series photons before they convert to Lya photons 
are negligible, because of the very large number of 
scatterings that the Lya photons undergo). This 
intensity is given by 

J r z e dF(z') \l + z' 1 . , , 

-=- = E a v a / E — — v a /E a dz , 

Jo Jz dz' [ 1 + z 

(31) 
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Fig. 12. — Ratio J/ Jo (lower curves) and y a (up- 
per curves) as a function of redshift. Solid curves 
are for the Press-Schechter mass function, dashed 
curves are for the Shcth-Tormen mass function. 



where E a = E(v a ). If the spectrum of first stars 
E(v) is approximated as flat near the Lya fre- 
quency (which is a good approximation for mas- 
sive, metal-free stars, which have surface temper- 
ature of 10 5 K), then the expression is simplified 
to 

i- = E a u a [F{z) - F{zp)] . (32) 
Jo 

We obtain an estimate for the value of E a as fol- 
lows: a typical halo with baryonic mass of 10 5 M Q 
forming one star with 1O 2 M uses a fraction 1CP 3 
of its mass to form stars. If ~ 50% of this mass is 
hydrogen that is burned to helium in the stellar in- 
terior during the star lifetime, the energy released 
per baryon in the halo is e& = 7MeVx 10~ 3 x 0.5 = 
3.5 keV. Assuming a blackbody spectrum with 
photosferic temperature T, we find 
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1 



hv a 

7T 4 exp(hu a /kT) - 1 \~~kT 



(33) 

For T — 10 5 K, this expression gives E a v a ~ 45. 
We use this value to compute the background in- 
tensity. 

The fraction of mass bound in halos above the 
minimum mass to form stars i s computed for the 
two halo models we use, i.e., the[P ress fc Schechterl 



F = erfc(v / s/2) , 
and the lSheth fc Tormenl l| 19991 ) model, 
0.4 



F = 0.4 



1 + 



,0.2 



(34) 



(35) 



where s = S^(z)/a 2 (M min ). The resulting Lya 
intensity is plotted as a function of z in Figure [T2l 
The corresponding y a , which gives the strengh of 
the spin-kinetic coupling by Lya scattering, is also 
shown in the same figure. 

The figure shows the Lya background intensity 
becomes important at very high redshift. For the 
Sheth-Tormen model, y a > 1 at z < 28, while 
for the Press-Schechter model, this is only slightly 
delayed (at z = 25). 

It may seem surprising that the background be- 
comes important at such an early stage in the for- 
mation of the first generation of stars, because as 
mentioned earlier the flux of local Lya photons in 
a Lya sphere is actually dominated by the sec- 
ondary Lya photons produced from the ionization 
by X-rays, whereas the Lya background is con- 
tributed only by the direct ultraviolet emission 
from the distant stars. But this is compensated for 
by the large number of first stars that are already 
visible from any location out to a large distance in 
the universe at these redshifts. This can be seen 
considering, as an illustration, that each star in 
the universe emits a luminosity Li in X-rays that 
are absorbed at radius ~ r a and have their en- 
ergy converted into "injected" Lya photons, and 
a luminosity L c of "continuum" photons in each 
frequency range Av jv a = Hr a /c, which enter the 
resonance Lya line as they propagate through a 
distance equal to the size of the Lya sphere. Then, 
the local flux of Lya photons from the star (dom- 
inated by the injected ones) is Fi = Li/ {Aitr 2 ^), 
and the background flux (determined by the con- 
tinuum photons) is F^ — n a L c r m , where n a is the 
number density of Lya spheres with first stars of 
similar luminosity, and r m is the maximum dis- 
tance out to which these first stars can be seen. 
In practice, r m is determined by the extremely 
fast evolution of the density of first stars (which 
parallels that of the background intensity shown 
in Fig. [T2|) . which doubles every Az ~ 2. Hence, 
r m ~ cH^ 1 Az I '(1 + z). Defining the filling factor 
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(36) 



of the Lya spheres / = 47rn Q r^/3, we find that 
F\ _ lj r a 1 
F b L c r m 3/ 

For the first stars at z ~ 30 and at r a — 10 kpc, 
Li/L c ~ 100 (as determined from a blackbody 
spectrum with T ~ 10 5 K), and r m /r a ~ 300, 
so the background of continuum photons becomes 
equal to the locally produced injected photons 
when the filling factor of the Lya spheres is about 
10%. The fluctuation in Lya photon flux as a 
function of smoo t hing sc ale was also discussed by 
Barkana fc Loebl (|2005af) . 

Hence, the background starts diminishing the 
contrast of the Lya spheres when their filling fac- 
tor becomes large. Of course, when the filling 
factor becomes large the background of injected 
Lya photons generated by X-rays from other stars 
(which we have not included here) should also be- 
come important. The importance of other stars 
in the neighborhood is greatly increased by the 
large biasing of the first stars. As we can see in 
Fig. [Til we typically expect to find several other 
halos capable of forming metal-free stars within 
the distance of the physical size of the Lya sphere 
(~ 10 kpc) at z < 25. This means that in practice, 
at z < 25 some of the Lya spheres may actually be 
produced by "clusters" of first stars, where several 
stars contribute to the Lya fluxes of large blobs 
yielding 21-cm absorption. 

4.2. The X-ray background 

In addition to the soft X-rays from the stellar 
photosphere, the intergalactic medium may also 
be heated up by a hard X-ray background, pro- 
duced for example by a population of early X-ray 
binaries after some of the metal-free stars have col- 
lapsed into black holes, or by supernova explosions 
and the relativistic electrons generated in super- 
nova remnants. The hard X-rays heat the atomic 
medium globally, because their mean free path is 
large compared to the typical separation between 
Lya spheres (as opposed to the soft X-rays from 
the metal-free main-sequence stars, which do not 
heat the medium substantially beyond the radius 
of a Lya sphere as discussed in §2). For local star 
bursts, it has been esti mated th at the energy re- 
lease in hard X-rays is (Oh 200 ll ) 



or about ex ~ 1 keV per baryon that forms stars. 
At high redshift the X-ray emission may increase 
because of the higher temperature of the CMB 
whic h provides seed photons for Compton scatter- 
ing dCerJl2003h . 

In practice, only a small fraction of the energy 
emitted in hard X-rays will be used to heat the 
IGM. At very high frequency, the universe is trans- 
parent to X-ray photons and so most of the emit- 
ted energy will be rcdshifted rather than absorbed 
by an atom. At redshift z = 30, the universe is 
opaque below a photon energy of ~ 2 keV. Only a 
small fraction of the energy of the absorbed pho- 
tons will be converted to heat (the rest being used 
for collisional ionizations and excitations) . If r\ e f f 
is the fraction of the emitted X-ray energy con- 
verted to heat, and /* is the fraction of baryons 
in collapsed halos with M > M m i n (z) that form 
stars, then the temperature evolution of the gas is 
described by (neglecting other heating and cool- 
ing mechanisms such as shock heating, Compton 
effect, etc.) 



dT _ F f*v eff e x | 2T 
dz fiks 1 + z 



(38) 



Lx 



5 x 10 40 erg 



SFR 

M Q yr- 



(37) 



To see a characteristic example for the gas tem- 
perature evolution, we integrate this equation us- 
ing /*?7e/ / = 2 x 10 -5 (a reasonable value obtained 
for /, ~ 2 x 10~ 3 for a typical metal-free star 
formed in a halo with ~ 10 5 M Q of baryons, and 
a fraction ry e // ~ 10~ 2 if 10% of the hard X-ray 
energy is emitted in the relevant range 0.2 — 2 keV 
and 10% of this energy is converted into heat). 
We start at z = 60 with a mean gas tempera- 
ture of 66.6 K as given by th e code RECFAST 
( Seager. Sasselov k, Scott 2000). The results are 
shown in Fig. [13] The CMB temperature and 
the adiabatically evolved gas temperature are also 
shown. The evolution of the mean temperature 
of the IGM depends sensitively on the hard X-ray 
emission. For ex ~ 1 keV per baryon, the mean 
gas kinetic temperature is below that of the CMB 
down to z ~ 11, but for e x ~ lOkeV, the gas is 
heated above the CMB temperature at z ~ 15 — 17 
in our model. 

We plot the profiles of the Lya sphere in the 
presence of hard X-ray heating in Fig. 1 141 with 
ex = lOkeV per baryon, f*r] e ff = 2 x 10~ 5 , and 
using PS halo model. At high redshift (z = 30 
to 25), when X-ray heating is still not significant, 
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Fig. 13. — Temperature evolution of the IGM. 
Solid curves show Tqmb and the adiabatically 
evolved gas temperature. Short and long-dash 
curves show the IGM heated by X-rays, with the 
PS (short dash) and ST (long dash) model calcula- 
tions. Two sets of curves are plotted, for e = 1 keV 
and lOkeV respectively. 



the profiles are similar to those shown in Fig. [4] 
At lower redshifts, however, the gas is gradually 
heated and therefore the 21-cm absorption be- 
comes gradually weaker. For e = 10 keV, the Lya 
sphere would of course appear in emission every- 
where for z < 20 because the gas kinetic temper- 
ature is globally heated above the CMB temper- 
ature (see the right hand panel of Fig [TTJ which 
shows a 200 Mq star at z = 15 where the mean 
temperature of the gas is already heated above the 
CMB temperature). 

It needs to be born in mind that the history 
of heating of the IGM is completely dependent on 
the type of X-ray sources that appeared among 
the first population of stars. These are likely to 
be X-ray binaries made by the first stellar black 
holes, in stars that are formed as close binaries. It 
is also possible that a population of miniquasars 
made when gas can cool in the vicinity of these 
first black holes emitted X-rays efficiently. The pa- 
rameters we have chosen here for the X-ray emis- 
sion should be considered as no more than an il- 



Fig. 14. — Brightness and temperature profiles of 
Lya spheres, as in Fig.0]except with X-ray heating 
of 10 keV per baryon. 



lustrating example, because we cannot predict the 
abundance of these first X-ray sources. The possi- 
bility of observing the transition from Lya spheres 
with a strong absorption signal to a weaker emis- 
sion signal at lower redshift would provide a very 
useful tool to study the thermal evolution of the 
atomic IGM and the nature of the first sources in 
the universe. 

Our simple estimate of the Lya and X-ray 
background are in gene ral agreemen t with other 



analyt ical models, e.g. [Se thi (2005J); iFurlanetto 



(|2006l ). However, ISethil (|2005l ) noted that the 



the actual process of X-ray heating may be an 
inhomogeneous one. Some part of the IGM gas 
would be heated earlier and to higher tempera- 
ture than other parts. With such temperature 
variations, one would observe that for the same 
redshift, some Lya spheres have only weak emis- 
sion signal, while others have strong absorption 
signal. Observation of such phenomena would re- 
veal tlie_tenipj3rja : ^^ the IGM (see 
also IPritchard & Furlanettoll2007[) . 
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5. Observability 

So far we have presented a theoretical investi- 
gation of the properties of what we have named " 
Lya spheres" , the first regions of the universe illu- 
minated by Lya photons from a metal-free star or 
a "cluster" of such stars (where the clusters would 
be caused by biasing in the cosmological distribu- 
tion of star- forming halos) . Would it be possible to 
observe this object directly with the incoming 21- 
cm instruments, such as 21CMA/PAST, LOFAR, 
MWA, and SKA? Or, if not these instruments, can 
we design one which could do it? What would be 
the characteristics of an instrument optimized for 
detecting the Lya sphere of the first stars? 

If the instrument has a response function 
R(v,9), with / dvd8 2 R(v, 6) = 1, then for a pixel 
centered on the Lya sphere, the resulting signal is 



T 



dvd8 2 R(v,9)ST(v, 9) 



(39) 



We can model R(v, 9) as a direct product of a top- 
hat spectral function with a bandwidth Av and a 
Gaussian beam of half-beam width O. Then, 



ST = 



1 



27re 2 A;X (2) _A, 



dv I d 2 9 5T(v,( 

IQ 



,-6 2 /2& 2 



(40) 



For an interfero meter array, the real space pixel 
noise is given by (iThompsom Moran. fc Swenson 



20011: ISantos. Coorav fc Knoxll2004l ) 



AT 



IJvt 



(41) 



where T sys is the system temperature, 5v is band- 
width, t the observation time, and 



fcov A dishAdish / ' -^-total 



(42) 



is the covering factor of the array. Here, Adish is 
the area of an antenna, and Atotai is the total area 
of the square made by the base line of the array. 
The angular resolution of the array is 



e 



-4, 



d i s h 



A 



total 



1/2 



A 
L 



(43) 



where L is the baseline, and A the redshifted wave- 
length. Thus, for a fixed value of NdishAdish^ the 




Bandwidth (MHz) 



Fig. 15. — The absolute value of the signal \6T\ 
as a function of beamwidth Q and bandwidth Av, 
for a star at z = 20, with a mass of 400 Mq and 
an age of 1.5 Myr. 



covering factor decreases with resolution. The sig- 
nal to noise ratio of the measured radio brightness 
within a resolution element of the array is 



SNR = VAvif c , 



ST 
'T~ 



(44) 



The system temperature is dominated by the 
galactic synchrotron background, which scales as 
2,-2.5 Qr j n t erms f rcdshift, 



1 + z 



2.5 



T sys ~ 2000K ■ (45) 

The signal to noise ratio at redshift z is then 
-2.5 / \ \ 1/2 



SNR ~ 3/ c . 



i±fV ( Ai/,t V f ST 

21 J VlOkHz-yrJ V 10mK 

(46) 

At redshift z, and for a flat space geometry, a 
comoving distance dr is related to angular size d9 
and frequency variation dv by 



d9 = dr I r, dv = vq 



Hdr 

c(l + z) 2 ' 



(47) 



where r — Jcdz/H. The comoving distance is 
related to the physical distance by r = (1 + z)r p . 
For the Lya spheres, the relevant physical scale is 
~ lOkpc. For the set of cosmological parameters 
adopted by us, a physical distance of 10 kpc at 
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Fig. 16. — The signal to noise ratio SNR as a 
function of beamwidth and bandwidth Ai>, for 
a star at z = 20, with a mass of 400 Mq and an 
age of 1.5 Myr, assuming f cov — 1 and t = lyr. 



redshift z = 20, 30 corresponds to an angular scale 
of 3.9, 5.3 arcsec and a frequency interval of 7.6, 
9.2 kHz, respectively. 

In a resolved observation, the maximum ab- 
sorption signal strength can reach a value close 
to ST — -200mK (Fig. HJ. However, in order to 
resolve the Lya sphere, a beam of < 4" and 
frequency bandwidth v < 8 kHz are required. For 
observations of the redshifted 21-cm line, the base- 
line for achieving an angular resolution 9 is given 

by 

+ (48) 



910 



21 



r 



km 



For a beamwidth and bandwidth similar to the size 
of the Lya sphere, one obtains the best signal-to- 
noise ratio for a single brightness measurement. 
Assuming that the center of the beam and fre- 
quency bin is centered on the star, we plot the 
observational signal \ST\ of the Lya sphere made 
by a 400M© star at z = 20 of age t = 1.5Myr, 
as a function of the beamwidth and bandwidth 
Av in Fig. 1151 For this star, the maximum emis- 
sion brightness occurs at 6.9 kpc away from the 
star, and the maximum absorption is located at 
15 kpc (ST = -190mK), but \ST\ > 5mK up to 
33 kpc. The signal decreases more rapidly with 
than Av, which is not unexpected given that the 
included volume is cx 2 Aza There are some wig- 
gles of ST in the Av direction, probably because 



as we vary the size of the frequency window, we 
are sampling different structure of the Lya sphere. 

We plot the corresponding signal to noise ratio 
(SNR; see eq. [55]) in FigureUHl assuming f cov — 1 
and t = lyr. Again, we ignore the second peak at 
very large bandwidth. To obtain a high confidence 
detection with SNR > 5(10, 20) under this circum- 
stance, one could choose an optimal bandwidth of 
30 kHz, and the optimal beamwidth (to obtain the 
highest SNR) is 20 (14, 10) arcsec, corresponding 
to a baseline of 45 (65, 91) km . The area corre- 
sponding to this baseline would have to be covered 
with dipole receivers in order to achieve f cov ~ 1. 
The requirements become more stringent with in- 
creasing redshift, owing to the increasing bright- 
ness of the synchrotron foreground, and as we have 
seen in §4 it may be necessary to go up to z ~ 30 
to see a pristine Lya sphere with a signal that is 
not yet weakened by a global Lya background and 
an IGM heated by hard X-rays. 

6. Conclusion 

We have investigated the 21-cm signature of the 
first stars. The coupling of the spin and kinetic 
temperature induced by the Lya photons gener- 
ated by a metal-free star should generate a "Lya 
sphere", a region that can be detected by a 21- 
cm absorption signal when the IGM has not been 
heated by X-rays prior to the formation of the star. 
The 21-cm absorption signal can reach a strength 
of 200 mK when the IGM temperature far from the 
star is at its minimum value reached by the adi- 
abatic cooling, and becomes weaker as the IGM 
is heated. Most of the Lya photons that couple 
the spin and kinetic temperatures are generated by 
collisional excitations of the high-energy electrons 
that are produced by soft X-rays emitted from the 
photospheres of the hot metal- free stars, rather 
than the direct photospheric emission of the star. 
These soft X-rays are not able to travel very far 
from the star and hence tend to heat the IGM only 
locally around each metal-free star that forms, but 
the IGM may be heated globally by other sources 
of hard X-rays, such as X-ray binaries formed by 
the black holes made by the first stars. 

The temperature and 21-cm absorption profile 
of a Lya sphere can be calculated as a function of 
the spectrum emitted by the star and its age in 
the idealized case of a uniform medium in Hub- 
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ble expansion around the star. However, owing 
to the highly biased distribution of the first stars, 
the region where a first star has formed may often 
be highly irregular, containing non-linear struc- 
tures on scales larger than the halo hosting the 
star. The first stars should also be highly clus- 
tered because of their biased distribution, so the 
high-redshift regions producing 21-cm absorption 
may often be produced by several first stars that 
are illuminating and heating a larger region with 
their X-rays and Lya photons. This may make 
these regions larger and therefore easier to detect. 

The detection of these earliest structures in the 
universe formed by the first stars is still very chal- 
lenging, mainly because of the high redshift at 
which they are likely to be present and the intense 
synchrotron foreground that must be subtracted 
to measure any CMB fluctuations. Filled aperture 
of at least ~ 50 km will be required for detection 
[| However, detecting these regions of Lya ab- 
sorption will represent a milestone for cosmology: 
they will constitute a clear indication that the first 
stars formed in the universe have been found, and 
allow us to measure many of their properties such 
as their mass distribution and spatial correlation. 
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Fig. 17. — A cross section map for a star of 
200 M Q and age 1.5 Myr at redshift 20(left) and 
15 (right). The box size is 40 kpc across (physical 
distance), and unit of temperature scale is mk. 
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